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Abstract 

Laser  cladding  is  an  additive  manufacturing  process  that  a  laser  generates  a  melt-pool  on  the 
substrate  material  while  a  second  material,  as  a  powder  or  a  wire  form,  is  injected  into  that  melt- 
pool.  Among  all  laser  manufacture  processes,  laser  cladding  offers  the  most  extensive  variety  of 
possibilities  to  alter  a  component  at  its  surface.  Despite  immense  potentials  and  advancements, 
the  process  model  of  microstructure  evolution  and  its  coupling  with  macro  parameter  of  laser 
cladding  process  has  not  been  fully  developed.  To  address  this  issue,  a  process  model  of 
microstructure  evolution  has  been  studied  by  utilizing  a  phase-field  method.  Phase-field  method 
has  become  a  widely  used  computation  tool  for  the  modeling  of  microstructure  evolution  with 
the  advantage  of  avoiding  of  tracking  the  interface  explicitly  and  satisfying  interfacial  boundary 
conditions.  In  present  work,  the  numerical  solutions  of  a  phase-field  model  are  analyzed.  The 
connection  of  macro-process  and  microstructure  evolution  is  examined  by  considering  the 
relationship  of  macro-  and  micro  parameters.  The  effects  of  thermal  noise  and  melt  undercooling 
on  the  final  microstructure  have  also  been  studied.  The  prediction  results  are  compared  with 
other  researchers’  results  and  good  agreement  was  found.  Different  solidification  morphologies 
of  different  locations  in  the  melt  pool  are  also  investigated.  It  was  found  that  not  the  mass 
transfer  but  the  heat  transfer  in  the  melt  pool  dominates  the  solidification  process. 
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I.  Introduction 


Laser  cladding  is  an  additive  manufacturing  process  that  a  laser  generates  a  melt-pool  on  the 
substrate  material  while  a  second  material,  as  a  powder  or  a  wire  form,  is  injected  into  that  melt- 
pool.  The  cladding  alloy  can  form  a  strong  bond  to  the  substrate,  with  a  minimum  melting  of  the 
substrate.  Among  all  laser  processes,  laser  cladding  offers  the  most  extensive  variety  of 
possibilities  to  alter  a  component  at  its  surface.  For  example  it  can  be  used  for  protection  of 
materials  against  the  wear,  corrosion  and  oxidation,  and  for  the  refurbishing  a  high  cost  industrial 
product  by  using  a  coating  with  improved  properties. 

During  the  laser  cladding  process,  several  complex  phenomena  occur  as  the  laser  interacts  with 
both  incoming  material  and  the  substrate  material.  These  phenomena  include  thermal  transport, 
fluid  flow,  mass  transport  and  solidification,  and  so  on.  Heat  and  fluid  flow  have  been  studied 
numerically  and  experimentally  by  other  researchers  [1,2].  There  are  many  researchers  who 
observed  the  microstructure  evolution  of  the  laser  cladding  process  by  using  optical,  SEM,  TEM 
technique  [25-27],  however,  the  modeling  of  microstructure  evolution  of  the  solidification 
process  during  the  laser  cladding  has  not  been  fully  developed  and  studied. 

Phase-field  method  has  recently  become  a  widely  used  technique  to  model  the  formation  of 
complex  interfacial  patterns  in  solidification  process.  It  has  been  the  main  topic  of  recent  review 
literatures  [4-5,  28].  With  this  model,  one  can  significantly  enhance  our  theoretical 
understanding  of  the  solidification  process,  especially  the  microstructure  evolution  in  the  metals, 
including  pure  substances  and  alloys.  In  the  phase-field  model,  the  interface  between  the  solid 
and  liquid  is  considered  as  a  diffuse  interface  instead  of  a  sharp  interface.  The  phase  field 
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variable  (j)  is  introduced  to  indicate  the  state  of  the  phase.  The  quantity  (j)  is  a  continuous  variable 
that  takes  a  constant  value  in  the  bulk  phases,  i.e.,  -1  for  liquid  and  1  for  solid,  and  changes  from 
-1  to  1  steeply  through  the  interface.  The  evolution  of  (])  is  governed  by  a  phase-field  equation, 
which  is  coupled  with  thermal  diffusion  equation  and  composition  equation  in  the  phase-field 
model  to  describe  the  transport  phenomena  in  the  solidification  process.  The  most  valuable 
property  of  the  phase-field  model  is  the  fact  that  there  is  no  need  to  explicitly  track  the  interface 
or  even  provide  interfacial  boundary  conditions  [6].  So  it  is  easy  to  implement  numerically 
compared  with  other  interface  front  tracking  method. 

Wheeler,  Boettinger  and  McFadden  developed  the  so-called  WBM  phase-field  model  to  deal 
with  alloy  solidification  [7,8].  This  model  was  derived  in  a  thermodynamically  consistent  way 
and  used  to  study  the  solute  trapping  of  the  isothermal  binary  alloy  in  rapid  solidification  [9].  It 
was  shown  that  the  solute  trapping  occurs  when  the  solute  diffusion  length  is  comparable  to  the 
diffuse  interface  thickness.  This  model  has  been  widely  used  in  many  applications,  including 
dendrite  growth,  solute  trapping,  eutectic  solidification,  sidearm  branching,  etc  [4].  Kim  et  al. 
[10]  proposed  another  phase-field  model  by  using  a  different  definition  of  the  free  energy.  This 
model  was  used  to  study  the  binary  alloy  solidification  and  dendrite  growth.  Their  results  are  in 
good  agreement  with  the  experimental  observations.  Wheeler  et  al.  [7]  have  shown  that  the 
interface  thickness  must  be  smaller  than  the  capillary  length  for  the  solution  to  converge  to  the 
sharp-interface  limit.  Karma  and  Rappel  [11]  reconsidered  this  issue  and  proposed  a  thin 
interface  analysis.  They  found  that  the  interface  thickness  only  needs  to  be  small  compared  to  the 
“mesoscale”  of  the  heat  and/or  solute  diffusion  field.  Their  asymptotic  analysis  greatly  enhanced 
the  computational  efficiency  of  the  phase-field  simulation  and  they  made  the  first  time  fully 
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resolved  eomputations  for  three-dimensional  dendritie  growth  [12,  13].  Beekermann  et  al.  [6,  14, 
15]  further  eondueted  mueh  effeetive  work  in  this  area.  They  presented  a  simpler  derivation  of 
the  phase-field  equation  starting  from  the  elassieal  veloeity-dependent  Gibbs-Thomson  interfaee 
eondition.  They  also  studied  the  pure  metal  dendritie  growth,  binary  alloy  solidifieation  eoupled 
with  heat  and  solute  diffusion,  and  the  solidifieation  phenomena  with  fluid  flow  [15]. 

Very  few  models  for  mierostrueture  evolution  of  laser  eladding  proeess  have  been  reported  [3, 
24],  and  almost  all  the  above  models  [6-15]  were  used  to  deal  with  the  ideal  proeess  (pure  metal, 
isothermal  binary  alloy,  ete).  In  this  paper  a  two-dimensional  phase-field  model  based  on  Karma 
and  Beekermann’ s  analysis  [13,  14]  was  developed  to  explore  the  mierostrueture  evolution  in  a 
real  physieal  proeess,  laser  eladding.  Anisotropy  effeet  was  ineluded  in  this  model,  while  the 
interfaee  kineties  was  negleeted  by  properly  ehoosing  the  ealeulation  parameters.  This  phase- 
field  model  was  eoupled  with  the  maero-proeess  model  and  served  as  miero-seale  model  to 
simulate  the  laser  eladding  proeess. 

This  paper  is  struetured  as  following;  The  maero-seale  model  of  laser  eladding  proeess  is 
presented  firstly.  This  model  is  used  to  get  the  maero-temperature  field  and  other  related 
physieal  parameters,  and  followed  by  the  phase-field  model.  The  linking  between  the  maero- 
seale  model  and  phase-field  model  is  earried  out,  eonsidering  the  relationship  of  the  maero-  and 
miero-seale  proeess  parameters.  The  phase-field  model  is  validated,  eomparing  the  interfaee 
shape  and  dendrite  tip  veloeity  with  other  researehers’  results  for  pure  metal  dendrite  growth. 
Next,  a  typieal  two-dimensional  binary  alloy  dendrite  growth  is  presented  for  the  ease  of  eoupled 
heat  and  solute  diffusion.  Then,  the  effeets  of  thermal  noise  and  melt  undereooling  on  the  final 
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microstructure  have  been  studied.  Different  solidification  morphologies  of  different  loeations  in 
the  melt  pool  are  also  investigated. 


II.  Macro-scale  Model 

The  maero-scale  model  was  used  in  this  study  to  get  the  maero-seale  temperature  field,  whieh 
will  be  used  in  the  phase-field  ealeulation.  The  transport  phenomena  of  melt  pool  were  studied, 
ineluding  the  melt  pool  shape  and  size,  and  so  on.  Most  work  of  this  part  is  based  on  the 
previous  study  of  the  author’s  researeh  group  [16]. 

The  basie  assumptions  of  this  model  are  as  follows: 

1 .  A  laser  beam  having  a  defined  power  distribution  strikes  the  surfaee  of  an  opaque  material  of 
finite  depth,  and  length. 

2.  Only  part  of  the  energy  is  absorbed  by  the  workpiece  and  powder.  Studies  show  that  the 
amount  of  laser  energy  absorbed  by  the  different  materials  is  37%-60%  [17,  23]. 

3.  The  absorbed  energy  induces  surfaee  tension  driven  flow  due  to  the  high-temperature  gradient. 

4.  The  liquid  metal  is  Newtonian,  so  the  Navier-Stokes  equation  is  applieable. 

5.  All  properties  of  the  liquid  and  solid  metal  are  eonstant,  independent  of  temperature.  (This 
allows  simplifieations  of  the  model;  however,  variable  properties  ean  be  treated  with  slight 
modifieations.) 

A.  Governing  Equations 

The  eontinuity  equation,  momentum  equation  and  energy  equations  based  on  eontinuum 
formulation  are  modified  and  used  in  the  present  study.  They  are  given  below: 
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Continuity  equation; 


-(p)  +  V-(pV)  =  0  (1) 

dt 


Momentum  equation: 


d  p  dp  p 

—  (Pm)  +  V-(pVm)  =  V-(p,— Vm)-  — -  — —  (m-m.)  (2) 

dt  p^  dx  K  Pi 


d  p  dp  p 

-(pv)  +  V-(pVv)  =  V-(p,^Vv)-^-^^(v-v^)  +  p.g  (3) 

dt  ^  Pj  dy  K  Pi 


Energy  equation: 


d(ph)  , 

+  V  •  (pYh)  =  V  •  (k  YT)  -  V  •  ipih,  -  h)(Y  -  VJ)  (4) 

dt  ‘ 


Evolution  of  free  surfaee: 


—  +  (Y-Y)F  =  0  (5) 

dt 

Note  that  the  momentum  and  energy  equations  eontain  the  phase  interaetion  terms.  In  the 
equations  (l)-(4),  the  subseript  1  and  s  denote  the  liquid  and  solid,  respeetively.  The  density  p, 
thermal  eonduetivity  kt  are  defined  as  the  average  value  of  liquid  and  solid  state  based  on  the 
solid  volume  fraetion  gs  and  liquid  volume  fraetion  gi.  The  speeifie  heat  ep,  veloeity  V  and  the 
enthalpy  h  are  defined  as  the  average  value  of  liquid  and  solid  state  based  on  the  solid  mass 
fraetion  fs  and  liquid  mass  fraetion  f.  u  and  v  are  the  X,  Y  eomponents  of  veloeity.  Us  and  Vs  are 
the  tangential  and  normal  veloeity  eomponents  at  the  free  surfaee.  K  is  the  permeability  of  the 
two-phase  mushy  region  whieh  is  modeled  as  a  porous  media,  pi  is  the  dynamie  viseosity. 


In  equation  (5),  E  expresses  the  fraetional  volume  oeeupied  by  fluid  in  that  eomputational  eell.  F 
takes  a  value  of  unity  at  a  eell  in  spaee  oeeupied  by  fluid  and  a  value  of  zero  if  fluid  does  not 
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occupy  that  cell,  and  cells  with  F  value  between  zero  and  unity  are  partially  filled  with  fiuid  and 
identified  as  surfaee  eells.  The  liquid-vapor  interfaee  is  loeated  at  the  F  =  0.5  eontour. 

B.  Boundary  Conditions 

With  the  liquid-vapor  free  surfaee,  the  shear  stress  and  normal  stress  balanee  is  eonsidered  to  get 
the  boundary  eondition  of  the  veloeity  at  the  free  surfaee.  The  energy  balanee  ineluded  the  effeet 
of  laser  heating,  and  heat  loss  from  eonveetion,  radiation  and  evaporation.  With  other  boundary 
surfaees,  the  veloeity  is  eonsidered  as  zero  and  energy  balanee  eonsidered  the  heat  loss  from 
eonveetion  and  radiation.  Details  of  the  formulations  for  the  boundary  eonditions  ean  be  found 
in  the  referenee  [16]. 


Ill,  Phase-field  Model 

One  maero-eell  adjaeent  to  solid-liquid  interfaee  will  be  ehosen  to  earry  out  the  mierostrueture 
evolution  ealeulation.  Thus  the  whole  ealeulation  domain  will  be  this  maero-eell.  The  behavior 
of  phase  field  variable  in  this  domain  is  governed  by  the  phase-field  equation,  whieh  is  eoupled 
with  equations  for  heat  and  solute  diffusion.  The  interfaee  between  liquid  and  solid  is 
represented  by  smooth  but  rapid  ehange  of  the  phase-field  variable.  Following  assumptions  were 
made  for  the  phase-field  model. 

1 .  Infiuenee  of  fiuid  flow  is  not  eonsidered  in  the  phase-field  model.  Only  a  very  tiny  maero-eell 
is  pieked  up  to  earry  out  the  miero-phase-field  ealeulation  sinee  the  effeet  of  fiuid  flow  is  not 
very  signifieant. 

2.  The  thermo-physieal  properties  sueh  as  thermal  eonduetivity,  speeifie  heat,  density,  and  so  on 
are  assumed  eonstant  and  equal  in  liquid  and  solid  state. 
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3.  Considering  the  solute  diffusivity  in  solid  state  is  very  small  eompared  with  that  in  liquid  state, 
the  solute  diffusion  in  solid  state  is  negleeted. 

4.  The  interfaee  kinetics  effect  is  neglected. 

5.  For  the  sake  of  simplicity,  the  partition  coefficient  of  the  solidification  process  is  assumed  to 
be  same  as  the  one  in  the  equilibrium  state. 


A.  Governing  Equations 

To  compare  the  result  with  other  researchers’,  the  dimensionless  temperature  and  composition 
are  employed  as  following: 

T  -T  -me 

d - ^  (6) 

Lie 

p 


X  = 


e-e^ 
(1  -  k)e^ 


(V) 


where  Tm  is  the  melting  temperature  of  pure  solvent,  m  is  the  liquidus  line  slope  of  the  alloy 
phase  diagram,  c„  is  the  value  of  composition,  c,  far  from  the  interface,  which  equals  the  initial 

composition  of  the  alloy,  L  is  latent  heat  of  melting,  k  is  the  equilibrium  partition  coefficient. 
Additionally,  the  following  dimensionless  variables  are  defined; 


and 


M  =  In 


2c /c^ 

1  +  A:  -  (1  -  k)(p 


(8) 


u  , 

e  - 1 
i-k 


(9) 


where  u  is  a  dimensionless  measure  of  the  deviation  of  the  chemieal  potential  from  its 
equilibrium  value  at  a  referenee  temperature,  and  U  is  used  to  ehange  the  variable  of  the 
eomposition  equation  and  considered  as  the  dimensionless  composition. 


The  phase-field  equation,  composition  equation  and  energy  equation  are  obtained  in  terms  of  the 
dimensionless  temperature  and  composition  as  follows: 


dt 


\  +  k  dU 
2  dt 


- {0  [l+(l-^)C/]}+  V  Z)^VC/  +  ^[1  +  (1-  k)m  —  ^ 

2  dt  I  2  2V2  dt  |V0| 


dO 

dt 


=  oV^d  + 


1  d(p 

2  dt 


where  x  is  the  relaxation  time  in  phase-field  model,  W  is  the  interface  thickness,  A,  is  a 


dimensionless  parameter  which  controls  the  strength  of  the  coupling  between  the  phase  and 
thermal/solutal  diffusion  fields,  and  D  is  the  solutal  diffusivity  in  the  liquid,  M  is  the  scaled 
liquidus  slope  which  is  defined  as; 


(10) 

(11) 

(12) 


M  = 


-m{l  -  k) 

Vs 


(13) 


To  eliminate  the  interface  kinetics,  x  is  chosen  to  be  a  function  of  U  in  the  phase-field  model, 
which  is  defined  by 

T  =To{— +  Mc„[l  +  (1-A:){/]}  (14) 

a 


where  Xo  is  the  measure  scale  of  time  in  current  study. 


B.  Anisotropy  Effect 
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To  include  anisotropy  effect,  we  choose  a  standard  fourfold  anisotropy  defined  by 


W(n)  =  W^a^(n)  t(«)  =  T„a/  (m) 

(15) 

V(/) 

(16) 

Combining  equations  (14)  -  (16),  we  ean  obtain  that 

T  («){—  + Mc„[l  +  (1-A:)[/]}  (17) 

a 

In  this  manner,  the  anisotropy  effeet  ean  be  ineorporated  and  the  interfaee  kineties  be  eliminated 
at  the  same  time. 


C,  Boundary  Conditions 

The  temperature  boundary  eondition  ean  be  found  from  the  maero-neighbor  eells  of  the  seleeted 
eell.  In  the  present  study,  the  Diriehlet  boundary  eondition  is  used  in  the  phase-field  caleulation. 
Beeause  it  is  assumed  there  is  no  fluid  flow  in  the  seleeted  eell  previously,  there  is  also  no  mass 
flow  between  the  whole  calculation  domain  and  neighbor  eells.  The  Diriehlet  boundary 
eondition  is  also  applied  for  eomposition  field. 

To  the  phase-field  variable,  a  zero-Newmann  eondition  is  employed  to  deseribe  the  boundary  of 
the  ealeulation  domain. 


IV,  Numerical  Method 

In  the  present  study,  the  maero-model  was  employed  to  ealeulate  the  temperature,  fluid  pressure 
and  veloeity  fields  of  the  whole  ealeulation  domain.  Each  macro-cell  has  four  state  variables: 
temperature,  pressure,  veloeity  and  fraetional  volume  (F).  One  maero-eell  near  the  solid-liquid 
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interface  is  chosen  to  carry  out  for  the  micro-scale  calculation  in  order  to  simulate  the  evolution 
of  solidification  microstructure.  This  macro-cell  is  re-meshed  making  much  finer  micro-cells. 
Each  micro-cell  has  three  state  variables:  temperature,  phase  field  and  concentration.  In  the 
temporal  domain,  there  are  two  time  steps:  At  for  macro-calculation  and  much  smaller  8t  for 
phase-field  calculation.  Each  At  is  again  divided  into  5  small  micro-time  steps  in  current 
calculation.  This  macro-  and  micro-scale  coupling  scheme  is  shown  in  Eig.  1. 

To  couple  the  two  different  scale  calculations,  the  temperature  of  the  neighbor  cells  of  the  micro¬ 
calculation  domain  is  applied  as  the  temperature  boundary  condition.  In  each  macro-time  step, 
this  boundary  temperature  is  considered  as  constant.  At  the  next  macro-time  step,  this  boundary 
temperature  will  be  changed  to  another  value,  which  is  determined  by  the  macro-scale  model. 

A  semi-implicit  SOEA-VOE  [18]  based  algorithm  was  used  for  the  solution  of  continuity  and 
momentum  equations  (l)-(3)  in  the  macro-scale  model.  The  thermal  energy  transport  equation 
(4)  was  solved  using  an  explicit  finite  volume  approximation.  A  fixed-grid  system  of  200  x  100 
uniform  grid  points  was  utilized  for  the  cladding  domain  of  4  mm  x  2  mm.  Each  grid  has  a  cell 
size  of  20  pm  x  20  pm.  An  explicit  finite  difference  scheme  was  used  to  solve  the  phase-field 
equation,  composition  equation,  and  energy  equation,  described  in  the  Equations  (10)-(12)  in 
micro-scale  model.  The  grid  numbers  used  in  this  study  are  either 1000  x  1000  or  1500  x  1500  in 
different  cases. 

It  is  assumed  that  the  nucleation  begins  at  a  certain  undercooling.  Initially,  the  micro-calculation 
domain  is  filled  with  an  amount  of  supercooled  melt.  Once  the  undercooling  reaches  the  given 
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value,  then,  the  nucleation  will  take  place  instantly  and  a  solid  seed  will  be  present  in  the  center 
of  the  calculation  domain,  followed  by  that  the  growth  process  of  the  solid  seed  will  be 
controlled  by  phase-field  model  and  appropriate  boundary  conditions. 

A,  Iteration  Procedures 

The  steps  of  iteration  procedure  for  the  macro-scale  calculation  are  as  following: 

(1)  Initialize  the  velocities,  pressure,  temperature  and  VOF  function. 

(2)  Solve  the  momentum  equations  utilizing  a  two-step  projection  method.  Surface  tensions  are 
solved  and  included  into  momentum  equations  as  body  forces.  In  this  step,  new  velocities  and 
pressure  are  solved. 

(3)  Advance  VOF  function  in  time  domain  to  update  new  free  surface. 

(4)  Solve  the  energy  equation  using  new  velocities  obtained  in  step  2. 

(5)  Check  the  undercooling.  If  it  reaches  the  given  value,  the  phase-field  calculation  will  begin. 
If  not,  go  to  step  6. 

(6)  Update  old  velocities,  pressure,  temperature  and  VOF  function  using  new  values,  and  back  to 
step  2  for  next  iteration  until  the  time  to  finish  calculation  is  reached. 

The  iteration  procedures  for  micro-calculation  are  as  following: 

(1)  Initialize  the  temperature,  composition  and  phase  field  variable. 

(2)  Solve  the  phase-field  equation  utilizing  an  explicit  finite  difference  method. 

(3)  Solve  the  composition  equation  using  new  phase-field  variables  obtained  in  step  2. 

(4)  Solve  the  energy  equation  using  new  phase-field  variables  obtained  in  step  2. 
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(5)  Update  old  temperature,  eoneentration  and  phase  field  variable  using  new  values,  and  baek  to 
step  2  for  next  iteration  until  the  phase-field  or  eomposition  field  reaeh  the  boundary  of 
ealeulation  domain. 

B,  Numerical  Stability 

With  the  maero-seale  ealeulation,  the  numerieal  differenee  equations  are  subjeeted  to  linear 
numerieal  stability  conditions  that  are  detailed  in  [19]. 

To  the  phase-field  model,  the  time  step  is 

dt  =  (18) 

where 

(19) 

2a  Wf, 

St,  is  derived  from  thermal  energy  equation  and  St^  from  phase-field  equation.  It  is  difficult  to 
determine  the  stable  time  step  for  composition  equation.  But  from  numerical  test,  it  is  found  that 
the  time  step  determined  by  equation  (18)  can  keep  the  composition  calculation  stable. 

V,  Simulation  Results  and  Discussion 

The  results  of  different  cases  are  presented  in  this  section.  First,  the  macro-scale  process  was 
simulated  using  the  macro-scale  model  developed  above.  To  validate  the  current  micro-scale 
model,  a  pure  metal  solidification  case  was  calculated  and  compared  with  Karma’s  result  [13]. 

A  binary  alloy  dendrite  growth  was  also  simulated.  Finally,  three  coupling  cases  were  calculated 
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with  different  loeations  in  the  melt  pool.  Simulations  of  all  the  eases  are  based  on  the  material 
properties  of  HI  3  tool  steel  and  proeess  parameters  shown  in  Table  1.  To  make  the  ealeulation 
simple,  the  material  used  in  this  study  was  assumed  as  a  binary  alloy. 

A,  Results  by  Macro-scale  Model 

To  carry  out  the  phase-field  calculation,  the  macro-scale  model  was  used  first  to  calculate  the 
whole  domain  temperature  field.  Several  snap  shots  of  the  calculation  results  are  shown  in  Fig.  2. 
Fig.  2(a)  shows  the  temperature  field  at  the  simulation  time  32  ms,  and  every  2  ms  until  38  ms, 
as  shown  in  Fig.  2(b)-(d).  Fig.  2  also  shows  the  powder  injection  process.  From  Fig.  2(a)  to  (d), 
it  can  be  seen  clearly  the  first  powder,  appeared  in  Fig.  2(a),  falls  into  the  melt  pool  as  shown  in 
Fig.  2(d).  It  is  also  observed  that  the  deeper  the  droplet  is  approaching  into  the  melt  pool  and 
pushing  the  fluid  away,  the  higher  the  speed  of  the  flow  is  at  the  top  surface. 

B,  Pure  Metal  Dendrite  Growth 

To  validate  the  current  phase-field  model,  a  pure  metal  solidification  case  was  calculated.  In 
pure  metal  case,  only  phase-field  equation  and  thermal  equation  need  to  be  solved.  And  the  Mc„ 
term  in  the  phase-field  equation  should  be  zero  in  this  pure  metal  solidification  case.  The 
dimensionless  undercooling  threshold  value  is  0.55  in  present  study.  According  to  Karma  [13], 
the  length  of  solidification  microstructure  is  scaled  by  Wo,  and  the  time  by  To.  Setting  Wo  =  1,  To 
=  1,  other  parameters  used  in  this  pure  metal  calculations  are  as  following:  X  =  6.3829,  a  =  4,  Ax 
=  Ay  =  0.4,  8t  =  0.008,  84  =  0.05.  The  grid  number  is  1000x1000.  Only  the  first  quarter  of 
calculation  domain  was  calculated  considering  the  symmetry.  For  the  comparison  with  Karma’s 
results,  the  interface  shape  is  calculated,  as  shown  in  Fig.  3,  and  the  interface  position  is  chosen 
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at  the  (j)  =  0  contour,  which  was  simulated  in  every  5000-iterations.  It  is  evident  that  this 
interface  shape  is  almost  identical  to  that  by  Karma’s  [13],  which  is  considered  that  this  model  is 
qualitatively  correct. 


The  dendrite  tip  velocity  was  also  calculated  in  this  pure  metal  case  according  to  the  following 
equation: 


d^/dt 

w 


(20) 


The  dimensionless  dendrite  tip  velocity  is  calculated  based  on  the  following  equation: 


Kip^O 

y  up  = - 


(21) 


where  do  is  the  capillary  length.  According  to  Karma’s  result  [13],  the  steady  state 
dimensionless  dendrite  tip  velocity  is  0.0170.  In  our  current  model,  the  converged  dendrite  tip 
velocity,  as  shown  in  Fig.  4,  is  0.0164,  which  is  pretty  close  to  Karma’s  result.  It  appears  that 
this  model  is  quantitatively  correct  and  the  difference  is  about  3.5%. 


C,  Binary  Alloy  Solidification 

In  order  to  compare  our  calculation  with  comparable  results  available  in  the  literature  [14],  it  is 
suitable  to  define  similar  parameters.  Thus,  the  Lewis  number  (Le),  which  is  defined  as  the  ratio 
of  thermal  diffusivity  a  to  solute  diffusivity  D,  is  chosen  and  used  for  this  comparison.  Based  on 
the  validation  of  pure  metal  solidification,  dendrite  growth  for  a  binary  alloy  case  was  simulated. 
The  simulation  parameters  are  chosen  as  following:  MC„  =  0.1,  k  =  0.15,  Le  =  20,  D  =  4.0,  8t  = 
0.0004.  Other  parameters  were  kept  as  same  as  pure  metal  dendrite  growth  case.  Fig.  5  shows 
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the  calculation  results  of  the  phase  field  ((])),  dimensionless  composition  field  (U),  and 
dimensionless  temperature  field  (0)  at  dimensionless  time  52  (i.e.  the  iteration  number,  as  8t  = 
0.0004  and  To=1,  is  52/0.0004=130000)  in  current  model.  The  final  solid/liquid  composition  can 
be  calculated  with  Equations  (8)  and  (9).  Comparing  Fig.  5(b)  with  5(c),  it  appears  that  the 
thermal  boundary  layer  thickness  is  much  larger  than  that  of  the  solutal  boundary  layer.  This 
arises  from  the  fact  that  the  thermal  diffusivity  is  20  times  larger  than  the  solutal  diffusivity.  The 
results  of  current  model  were  compared  with  Ramirez’s  [14].  Excellent  agreement  was  also 
found  except  the  range  of  composition  field.  It  was  explained  that  different  value  of  Eewis 
number  was  used  in  the  calculation,  which  is  Ee=20  in  current  model,  whilst  Ee=40  in  Ramirez’s 
calculation. 

D,  Results  for  Coupling  Cases 

Three  different  coupling  cases  based  on  the  location  in  the  melt  pool  were  investigated  in  this 
study.  The  cell  location  is  shown  in  Fig.  6.  All  other  simulation  parameters  were  kept  same  in 
these  three  cases. 

Case  1:  The  macro  cell  was  chosen  from  cell  location  1  in  Fig.  6,  which  is  located  on  the 
solid/liquid  interface  and  on  the  left  of  melt  pool.  Considering  the  workpiece  move  direction  (to 
the  left),  one  can  figure  out  that  there  will  be  certain  incoming  heat  transfer  from  right  to  left 
during  the  following  time  step.  This  can  be  shown  clearly  from  the  dimensionless  temperature 
evolution  of  Fig.  7.  Fig.  7  also  shows  the  microstructure  and  composition  evolution  of  cell 
location  1  from  dimensionless  time,  10  to  60.  An  interesting  observation  is  that  the  left  dendrite 
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tip  shows  a  greater  growth  veloeity  than  the  right  one  with  the  ineoming  heat  flux  from  right  to 
left  side. 

Case  2:  The  maero  eell  was  pieked  up  from  the  left  bottom  eorner  of  the  melt  pool,  whieh  is  eell 
loeation  2  in  Fig.  6.  The  microstrueture,  dimensionless  temperature  and  eomposition  evolution 
were  shown  in  Fig.  8.  It  ean  be  found  that  the  heat  flux  direetion  in  Fig.  8  is  along  the  diagonal 
of  the  ealeulation  domain,  from  the  right  top  eorner  to  the  left  bottom  corner,  and  that  the  left  and 
bottom  dendrite  tips  grow  faster  than  the  right  and  top  ones. 

Case  3:  The  macro  cell  was  selected  from  the  cell  location  3  shown  in  Fig.  6.  This  cell  locates  in 
the  right  front  of  the  melt  pool,  which  will  soon  be  heated  by  the  external  energy  source  (laser  in 
this  study).  Fig.  9  shows  the  corresponding  evolution  of  microstrueture,  dimensionless 
temperature  and  composition.  The  initial  seed  in  the  calculation  domain  undergoes  a  growth  in 
the  first  two  frames  and  subsequently  melt  in  the  remaining  frames.  The  interesting  point  is  that 
the  faster  growth  direction  in  the  first  two  frames  is  identical  to  the  heat  flux  direction,  which  is 
from  left  top  corner  to  the  right  bottom  corner. 

Based  on  the  observation  of  the  above  three  cases  which  share  same  calculation  parameters 
except  the  location  in  the  melt  pool,  one  can  reach  the  conclusion  that  it  was  the  heat  transfer,  not 
the  mass  transfer,  that  dominates  dendrite  growth  pattern  in  the  binary  alloy  solidification 
microstrueture  during  laser  cladding  process. 

E,  Effect  of  Thermal  Noise 


17 


It  was  reported  that  the  conserved  noise  in  heat  current  is  the  most  dominant  one  according  to 
Karma’s  study  on  dendritic  sidebranching  with  thermal  noise  [20].  In  the  present  study,  a 
different  method  [21]  is  employed  to  include  the  thermal  noise  effect.  To  simplify  the 
calculation,  only  the  pure  metal  case  was  calculated  in  this  section.  The  thermal  energy  equation 
with  thermal  noise  changes  to 

- aV  9  +  —  +  5’g(l  +  0)(1  -  0)  (22) 

dt  dt 

where  S  is  the  thermal  noise  strength,  ^  is  a  random  number  from  -1  to  1.  The  thermal  noise 
term  in  Equation  (22)  shows  that  the  thermal  noise  will  be  vanished  at  the  bulk  state  (solid  or 
liquid)  and  only  valid  at  the  diffuse  interface  region.  It  is  obvious  that  the  highest  thermal  noise 
is  present  if  the  (j)  equals  to  0,  where  the  solid/liquid  interface  is  defined  in  this  work. 

Fig.  10  shows  the  effect  of  thermal  noise  on  the  dendrite  growth  pattern,  varying  thermal  noise 
strength.  When  the  thermal  noise  strength  is  small  (S  =  4),  only  very  small  variation  of  the 
growth  pattern  can  be  observed  along  the  dendrite.  With  the  increase  of  the  strength,  the 
variation  becomes  larger  and  easier  to  observe.  It  can  be  observed  that  different  dendritic 
sidebranching,  showing  secondary  arm,  even  the  tertiary  arm,  appears  in  the  growth  pattern. 

Also,  it  is  observed  that  the  asymmetry  effects,  which  arise  from  the  random  number  (^)  in  the 
thermal  noise  term,  can  be  present  in  the  final  microstructure.  When  the  thermal  noise  strength 
is  much  higher,  i.e.,  S  =  25,  as  shown  in  Fig.  10  (e),  the  bifurcation  phenomenon  appear  in  one  of 
the  primary  dendrite  arm.  The  detailed  mechanism  of  this  phenomenon  is  beyond  of  this  study 
and  still  investigating,  thus  will  not  be  discussed  further. 
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F,  Effect  of  Undercooling 

The  effect  of  undercooling  on  the  dendrite  growth  pattern  of  the  solidification  microstructure 
was  also  investigated.  In  this  study,  all  other  parameters  were  kept  same  with  the  pure  metal 
case,  and  only  the  threshold  undercooling  was  varied  from  0.45  to  0.75. 

Fig.  1 1  shows  the  calculated  phase  field  with  different  undercooling.  It  is  clear  that  one  can  get 
sharp  and  smooth  dendrite  shape  if  we  choose  A  =  0.55.  If  the  undercooling  is  smaller  than  the 
value,  as  shown  in  Fig.l  1  (a),  the  primary  dendrite  arm  will  become  coarser.  Meanwhile,  if  the 
undercooling  is  greater  than  0.55,  as  shown  in  Fig.  1 1  (b),  there  will  be  some  bifurcation  in  the 
final  phase  field.  Increasing  the  undercooling  much  higher,  as  shown  in  Fig.  1 1  (c),  much  more 
bifurcation  will  be  appeared.  It  appears  that  the  dendrite  growth  pattern  varies  from  dendrite  to 
globular  as  the  undercooling  increases.  It  also  turned  out  that  the  CPU  time  for  the  simulation  is 
decreasing  dramatically  with  the  increase  of  the  undercooling.  Considering  the  same  domain 
size,  it  is  evident  that  the  dendrite  tip  velocity  increases  dramatically  with  the  increase  of  the 
undercooling.  This  phenomenon  was  also  reported  by  Zhao  et  al.  [22],  and  according  to  their 
analysis  the  dendrite  tip  velocity  is  proportional  to  in  2D  case  with  small  Peclet  number. 

VI.  Concluding  Remarks 

1 .  A  multi-scale  model  for  laser  cladding  process  has  been  developed  by  coupling  the  macro¬ 
scale  model  based  on  VOF  formulation  with  micro-scale  model  based  on  phase-field  model.  In 
phase-field  model,  the  effect  of  anisotropy  was  included  in  predicting  dendrite  growth  pattern  but 
the  interface  kinetics  was  neglected  to  simplify  the  modeling  work.  This  multi-scale  model  can 
be  used  to  effectively  simulate  solidification  microstructure  during  laser  cladding  process. 
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2.  To  verify  the  miero-scale  model  based  on  phase-field  approaeh,  pure  metal  dendrite  growth 
was  investigated.  A  good  agreement  was  found  with  the  results  reported  in  literature.  Dendritie 
growth  for  the  binary  alloy  was  also  simulated,  and  the  temperature  and  composition  distribution 
in  the  dendritic  growth  showed  similar  patterns  observed  by  other  researchers. 

3.  Phase-field,  dimensionless  temperature  and  composition  at  different  cell  locations  in  the  melt 
pool  during  the  laser  cladding  were  investigated  linking  macro-scale  model  with  micro-scale 
model.  It  was  found  that  the  heat  transfer,  not  mass  transfer,  dominates  dendrite  growth  pattern 
of  solidification  microstructure  for  the  binary  alloy  case. 

4.  The  effects  of  thermal  noise  and  undercooling  on  the  growth  pattern  of  solidification 
microstructure  were  also  studied.  It  was  found  that  some  secondary  arm  or  even  tertiary  arm  in 
the  dendrite  growth  process  could  be  side-branched  if  appropriate  thermal  noise  strength  was 
chosen.  Also,  if  the  undercooling  was  chosen  appropriately,  smooth  dendrite  pattern  with 
suitable  dendrite  tip  velocity  could  be  obtained. 
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anisotropy 

composition  of  material 

specific  heat  of  liquid 

specific  heat  of  solid 

composition  far  from  the  interface 

capillary  length 

solutal  diffusivity 

volume  of  fluid  function 

liquid  mass  fraction 

solid  mass  fraction 

liquid  volume  fraction 

solid  volume  fraction 

enthalpy 

liquid  enthalpy 

solid  enthalpy 

grid  number  in  x  direction 

grid  number  in  y  direction 

equilibrium  partition  coefficient 

thermal  conductivity 

permeability  of  the  two-phase  mushy  region 
latent  heat  of  melting 
Lewis  number 

latent  heat  of  fusion  of  material 

latent  heat  for  liquid-vapor  phase  change 

slope  of  liquidus  line 

scaled  liquidus  slope 

pressure 

thermal  noise  strength 

temperature 

liquidus  temperature 

melting  temperature  of  pure  solvent 

solidus  temperature 

time 

macro-time  step 
micro-time  step 

dimensionless  measure  of  composition 

velocity  in  X-direction 

velocity  in  Y-direction 

tangential  velocity  at  free  surface 

normal  velocity  at  free  surface 

velocity  vector 

dendrite  tip  velocity 


24 


Vtip 

w 

K 

A 

Ax 

Ay 

Subscript 

I 

s 

Greek  symbols 

a 

X 

£4 

4> 

Hi 

A 

e 

p 

Pi 

Ps 

r 

To 


dimensionless  dendrite  tip  veloeity 
interfaee  thiekness 
length  scale  in  eurrent  study 
undereooling 

mesh  spaeing  along  X  axis 
mesh  spaeing  along  Y  axis 

liquid  phase 
solid  phase 

thermal  diffusivity 
dimensionless  eomposition 
anisotropy  eoeffieient 
phase-field 

dynamie  viseosity  of  liquid  metal 

strength  of  eoupling  between  phase  and  diffusion  fields 

dimensionless  temperature 

density 

density  of  liquid 
density  of  solid 

relaxation  time  in  phase-field  model 
time  scale  in  eurrent  study 
Random  number  between  -1  and  1 
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Table  Captions 

Table  1.  Material  properties  and  process  parameters 
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Table  1.  Material  properties  and  proeess  parameters 


Property 

Symbol 

Value 

Powder  absorptivity 

a 

0.35 

Power  absorptivity 

p 

0.3 

Specific  heat  of  liquid 

Cpl 

800  J/(kgK) 

Specific  heat  of  solid 

Cps 

700  J/(kgK) 

Diameter  of  powder 

D 

100  pm 

Conductivity  of  liquid 

ki 

21  W/(mK) 

Conductivity  of  solid 

ks 

21  W/(mK) 

Latent  heat  of  melt 

Lm 

2.7x  lOM/kg 

Latent  heat  of  evaporation 

Lv 

6.34x  lOM/kg 

Powder  mass  rate 

nip 

5  g/min 

Power  of  laser  beam 

P laser 

1000  W 

Density  of  liquid 

Pi 

6900  kg/m" 

Density  of  solid 

Ps 

7800  kg/m" 

Radius  of  laser  beam 

R 

0.7  mm 

Liquidus  temperature 

Ti 

1700  K 

Solidus  temperature 

Ts 

1630  K 

Melting  temperature 

T™ 

1665  K 

Dynamic  Viscosity 

Pi 

6  X  10'"  kg/(ms) 

Velocity  of  laser  beam 

^  laser 

12.7  mm/s 

Powder  velocity 

V  powder 

0.2  m/s 
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Figure  Captions 

Fig.l  Scheme  of  maero/miero-seale  eoupling  method. 

Fig. 2  Temperature  field  of  whole  domain  (powder  veloeity,  0.2mA,  powder  diameter,  100/im, 
laser  beam  power,  lOOOW,  radius  of  laser  beam,  0.7  mm,  and  laser  absorptivity,  0.3). 

Fig. 3  Pure  metal  dendrite  growth:  Interfaee  shape  (eontour  (])  =  0)  shown  every  5000  iterations, 
simulation  parameters,  X  =  6.3829,  a  =  4,  Ax  =  Ay  =  0.4,  8t  =  0.008,  £4  =  0.05. 

Fig. 4  Convergenee  of  dimensionless  dendrite  tip  veloeity  as  inereasing  iteration  number  for  the 
simulation  of  pure  metal  dendrite  growth. 

Fig. 5  Caleulation  results  for  alloy  solidifieation  (at  dimensionless  time  of  52)  (MC„  =  0.1,  k  = 
0.15,  Le  =  20,  D  =  4.0,  8t  =  0.0004,  Ax  =  Ay  =  0.4,  £4  =  0.05)  (a)  Phase-field  distribution 
(b)  Dimensionless  composition  field  distribution  (e)  Dimensionless  temperature  field 
distribution. 

Fig. 6  Scheme  of  eell  loeation. 

Fig. 7  Mierostructure  evolution  of  eell  loeation  1  (MC„  =  0.1,  k  =  0.15,  Le  =  20,  D  =  4.0,  8t  = 
0.0004,  Ax  =  Ay  =  0.4,  84  =  0.05). 

Fig. 8  Mierostructure  evolution  of  cell  loeation  2  (MC„  =  0.1,  k  =  0.15,  Le  =  20,  D  =  4.0,  8t  = 
0.0004,  Ax  =  Ay  =  0.4,  84  =  0.05). 

Fig. 9  Mierostrueture  evolution  of  eell  loeation  3  (MC„  =  0.1,  k  =  0.15,  Le  =  20,  D  =  4.0,  8t  = 
0.0004,  Ax  =  Ay  =  0.4,  84  =  0.05). 

Fig.  10  Effeet  of  thermal  noise  strength  (S)  on  dendrite  growth  pattern:  (a)  S=4,  (b)  S=10,  (e) 
S=15,  (d)  S=20,  and  (e)  S=25. 

Fig.  1 1  Effeet  of  undereooling  (A)  on  dendrite  growth  pattern:  (a)  A  =  0.45,  (b)  A  =  0.65,  and  (c) 
A  =  0.75 
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Fig.l  Scheme  of  macro/micro-scale  coupling  method 
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(b)  t  =  34  ms 


(c)  t  =  36  ms 
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(d)  t  =  38  ms 

Fig.  2  Temperature  field  of  whole  domain 


(powder  velocity,  0.2m/s,  powder  diameter,  lOOjum,  laser  beam  power,  lOOOW,  radius  of  laser 


beam,  0.7  mm,  and  laser  absorptivity,  0.3). 
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Fig.  3  Pure  metal  dendrite  growth;  Interface  shape  (contour  (])  =  0)  shown  every  5000  iterations, 
simulation  parameters,  A.  =  6.3829,  a  =  4,  Ax  =  Ay  =  0.4,  8t  =  0.008,  £4  =  0.05. 


31 


Dimensionless  Tip  Velocity 


Fig.  4  Convergence  of  dimensionless  dendrite  tip  velocity  as  increasing  iteration  number  for  the 

simulation  of  pure  metal  dendrite  growth 
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Fig.  5  Calculation  results  for  alloy  solidification  (at  dimensionless  time  of  52) 
(MC„  =  0.1,  k  =  0.15,  Le  =  20,  D  =  4.0,  6t  =  0.0004,  Ax  =  Ay  =  0.4,  84  =  0.05) 

(a)  Phase-field  distribution,  (b)  Dimensionless  composition  field  distribution,  and  (c) 
Dimensionless  temperature  field  distribution 
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Fig. 6  Scheme  of  cell  location 
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Fig.  7  Microstructure  evolution  of  eell  loeation  1 
(MC  =  0.1,  k  =  0.15,  Le  =  20,  D  =  4.0,  8t  =  0.0004,  Ax  =  Ay  =  0.4,  S4  =  0.05) 
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Fig.  8  Microstructure  evolution  of  cell  location  2 
(MC„  =  0.1,  k  =  0.15,  Le  =  20,  D  =  4.0,  8t  =  0.0004,  Ax  =  Ay  =  0.4,  84  =  0.05) 
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Fig.  9  Microstructure  evolution  of  eell  loeation  3 
(MC„  =  0.1,  k  =  0.15,  Le  =  20,  D  =  4.0,  6t  =  0.0004,  Ax  =  Ay  =  0.4,  84  =  0.05) 
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(a)S=4 


(b)  S  =  10 


(c)  S=  15 


(d)  S  =  20 


(e)  S  =  25 


Fig.  10  Effect  of  thermal  noise  strength  (S)  on  dendrite  growth  pattern 
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Fig.  1 1  Effect  of  undercooling  (A)  on  dendrite  growth  pattern 
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